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ABSTRACT 



Iron spectral features are thought to be the best tracer of a progenitor of 



cn . gamma-ray bursts (GRBs). The detections of spectral features such as an iron 

f^^ . line and/or a Radiative Recombination edge and Continuum (RRC) were re- 

O , ported in four X-ray afterglows of GRBs. However their properties were different 

each other burst by burst. For example, Chandra observation of GRB 991216 
^ . reported both the strong H-like iron line together with its RRC. On the contrary, 

Yoshida et al. (2001) report only a detection of the strong RRC in GRB 970828 
with ASCA. Since it is difficult to produce the strong RRC, we have to consider 



O 

c/3 . special condition for the line and/or the RRC forming region. In this paper, we 

point out a possibility of a "non- equilibrium ionization state" for the line and the 



RRC forming region. 
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1. Introduction 

Gamma-Ray Bursts (GRBs) are remote at cosmological distances and tlius tlieir released 
energies in 7-ray photons are almost ~ 10^^ erg (Kulkarni et al. 1999). Although fireball 
models (Rees & Meszaros 1992; Piran 1998) or a cannonball model (Dar 2000) can explain 
many observational properties of GRBs, a progenitor of GRB is yet to be solved. An impor- 
tant key to probe a progenitor of GRB is the detection of iron spectral features in X-ray after- 
glows. Iron is the most abundant heavy element and fall in the middle of the observed energy 
range of X-ray satellites. Strong iron features can be produced in a dense gas environment. 
Therefore, the detection is the best evidence of a collapsing model of a massive star such as 
Hypernova/Collapsar or Supranova (Paczyhski 1998; Woosley 1993; Vietri et al. 1998). 

The iron features were already reported in four X-ray afterglows. The redshifted 
iron emission line was discovered in the X-ray afterglow of GRB 970508 with BeppoSAX 
(Piro et al. 1999), and the redshift was consistent with the distance of a host galaxy. An 
independent discovery of a redshifted iron emission line was reported for GRB 970828 with 
ASCA (Yoshida et al. 1999). Assuming a He-like Ka line, the spectral feature was once 
interpreted as an emission line with a redshift oi z = 0.33. However, the Keck observation 
of the host galaxy, which was discovered later in 1998, revealed the redshift oi z = 0.9578 
(Djorgovski et al. 2000). The discrepancy in distance with ASCA and moreover the tempo- 
ral detections of both forced us to doubt the reality of existence of the iron features before 
the confirmation by the Chandra detection. 

In 1999, Chandra observed GRB 991216 at 37 hours after the GRB with the ACIS- 
S/HETG instruments, and Piro et al. (2000) reported the detection of both the iron emission 
line together with the Radiative Recombination edge and Continuum (RRC) from fully 
ionized iron with high statistical significance of 4.7 a and 3.0 a respectively. The redshift 
of the line and the RRC were consistent with the host galaxy, and the line was broaden 
to (Tiine = 0.23 ± 0.07 keV suggesting a moving ejecta. Soon after the Chandra detection, 
Antonelli et al. (2000) also reported the probable iron emission line in the spectrum of 
GRB 000214 with ~ 3.2 a confidence level. This was the second detection of the iron with 
BeppoSAX. However the redshift of a host galaxy of GRB 000214 is unknown. Therefore, it 
remains a question whether the observed feature is the iron emission line or the RRC. 

Motivated by the RRC detection with Chandra, Yoshida et al. (2001 paper-I) apply the 
RRC for the ASCA data of GRB 970828 using the redshift oi z = 0.9578, and obtain the edge 
energy, which is consistent with 9.28 keV. This possibility was earlier suggested in the paper 
by Djorgovski et al. (2001). However, there was no iron Ka line detected at the expected 
energy from GRB 970828. The upper limit of the line flux is < 1.5 x 10^^ photons cm^^s^^ 
or 120 eV in equivalent width (EW). We must explain the no detection of Ka line. 
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Related to the iron features, we should note the fact that the iron features were found 
only in a small fraction (< 10 %) of X-ray afterglows, and sometimes the features were 
observed only during a certain interval. BeppoSAX has tried to search lines in eleven X-ray 
afterglow spectra, but found only one from GRB 970508 at the end of 1999 (Piro; private 
communication). Most X-ray afterglows did show only upper limits. Especially, Yonetoku 
et al. (2000) set an extremely low upper limit of almost 100 eV in EW for the bright 
GRB 990123 with ASCA. The published intensities of the observed iron Ka line and the 
RRC are summarized in table 1. There is a wide variety in the iron emission and the 
RRC intensity. In this paper, we try to explain these varieties of the iron features with 
the assumption that the line-emitting plasma state is in "non-equilibrium ionization (NET) 
state" which has a low electron temperature as compared to the ionization degree. 



2. Spectral Simulation for Non-equilibrium Plasma 

In paper-I, the observed flux of the RRC and the upper limit for the iron Ka line are 
given for GRB 970828 with 90 % confidence errors. So, we focus on the ratio of the integrated 
RRC flux to the iron intensity: -FRRc/-^iinc, which is free from the continuum. The observed 
ratio of: -FRRcZ-^une > 3.3 in 90 % statistical lower limit is very hard to reproduce. The case 
of a strong RRC without Ka line looks abnormal. In fact, BeppoSAX detected the strong 
Ka line and the Chandra observation showed both features of H-like iron with the flux ratio: 
-^RRcZ-^iine ~ 1- Therefore, we are forced to consider a different condition between ASCA 
and other results containing the strong Ka line. 

The no Ka line of highly ionized iron accompanied by the recombination edge gives 
a constraint on the possible emission mechanisms. If the line were produced due to ex- 
citation by electron-impact, the RRC would be hidden by a much more intense thermal 
bremsstrahlung. Also the synchrotron emission, which is thought to dominate the afterglow, 
would be significantly overlaid by the thermal bremsstrahlung. Therefore, we need a condi- 
tion that the iron is highly ionized but the involved electron energy is low. This is a radiative 
recombination in the NEI (Tg < T^) state or charge exchange must be responsible for the 
iron emission. 

In this section, we show spectra by numerical calculations in the NEI plasma state, not 
depending on the specific model. To explain the observed strong RRC without the Ka line 
of iron, we calculate the emissivity using the NEI plasma radiation code (Masai 1994). The 
code employs three mechanisms for the continuum, such as free-free emission, two-photon 
decay and radiative recombination. For line emissions, as well as excitation by electron- 
impact, fluorescence lines due to ionization and cascade lines due to recombination are taken 
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into account. Radiation properties of a plasma were described by two parameters of the 
electron temperature (Tg) and the ionization degree represented in units of temperature (T^), 
assuming a cosmical abundance. We study the emissivities in the range of 0.1 < kT^, < 10 keV 
and 0.1 < kTz < 100 keV in every 0.1 keV step. We show representative spectra in figure 1. 

The strong RRC compared with the Ka line can be formed only in the regime of T^ < Tz, 
recombining plasma condition. This condition is, however, attained by several situations as 
discussed later. We intend to find the condition of the plasma to account for the observed flux 
ratio (i^RRcZ-^iine), which is free from specific modelling with iron abundance, the emission 
measure, the geometry and so forth if the line and RRC are emitted from the same cite. 
Thus, for a given Tg and T^ a priori, we carried out calculations of the emissivity ratio in the 
above wide range of Tg - T^ space. 

Figure 1-c and 1-d are the simulated plasma emissivities with a cosmical abundance to 
explain the ratio: -FRRcZ-^iine of GRB 991216 and GRB 970828 respectively. The ratios of 
an integrated RRC flux to Ka lines of the simulated ionization state are -FRRc/-^ime ~ 1 and 
~ 4 for figure 1-c and 1-d respectively, and the ratios are within the observational values 
in 90 % statistical error. In the calculation of the ratios, the line components (mostly the 
blend of Fe and Ni) very close to the edge of the RRC are included to the RRC component, 
due to the limited energy resolving power of the SIS detectors onboard ASCA. 



3. The Reason for Strong RRC and Weak Ka Line 

To produce the strong RRC in quantum number n = 1, the iron must be almost fully 
ionized. The H-like Ka line, which was observed with Chandra, dominates other ionization 
states at a temperature of kTz > 20 keV. Above the temperature, Fe XXVII consists of 
more than 70 % of iron, thus we assume kTz > 20 keV in the following discussion. In a 
condition of high electron temperature of kT^ ~ kTz > 20 keV, i.e., equilibrium ionization, 
the capture rate of free electrons is small and the emissivity of the line and RRC also becomes 
small. Moreover, the free-free emission from high Tg electrons dominates at the hard X-ray 
band. Therefore, the RRC may not be observed because it can be obscured by the free-free 
component. 

The cross section of the electron capture into n th quantum state can be expressed as 

1 S kle 1 , _i , , 

(^nOC—{- + — ) (1) 

n^ 2 e n"^ 

where e is an ionization energy of 9.28 keV for H-like iron (e.g. Nakayama & Masai 2001). 
Thus, the best condition to form the strong RRC would be the case of kTe ~ e (<C kTz). In 
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such a plasma state, o"„ oc n~^ and then most of free electrons recombine directly into the 
ground state (n = 1), compared to the n > 2 levels. However, free-free emission dominates 
the continuum. 

With decreasing kTe, the recombination rate increases, while free-free emission becomes 
suppressed. Recombination into n > 2 increases relatively and produces line emission. Thus, 
the Ka line can be enhanced by cascades from n > 3 excited levels. This is the case for 
He-like Ka, but H-like Ka (Lya) is little affected; a considerable fraction comes to direct 
transition to the ground state. 

We summarize the above discussions in view point of the intensity of the RRC and Ka 
line. The plasma state with the strong RRC but the weak Ka line, which was observed with 
ASCA, is realized when kT^ is slightly less than e but in high kT^. The ratio of the RRC to 
Ka lines (-pRRcZ-^ime) from the numerical calculations is shown in figure 2 as the function 
of Te and T^. The peak of the ratio appears at around kT^ = 4 keV and kTz = 100 keV in 
the calculated range. The condition of -FRRcZ-^Une > 3.3 of ORB 970828 can be explained by 
this result. 



4. Discussion 

A high ionization degree of kTz ~ 100 keV and low electron temperature of kT^ ~ 1 keV 
are required when we reproduce the high ratio of the RRC to the iron line more than 3.3. 
This condition may be attained in the situations: (i) photoionizations by X-rays or (ii) rapid 
cooling due to rarefaction. 

In the case (i), also a fluorescence Ka line of energies 6.4 - 6.5 keV of partially ionized 
iron is likely accompanied. Especially, the photoionizations of a neutral circumstellar gas by 
the initial bright flash of GRB, iron lines in the low ionization states are expected but not 
observed. The line observed with Chandra was purely H-like (Piro et al. 2000). Therefore, 
the iron atoms should be fully ionized by the time of the observed iron emission. However, if 
the line and the RRC emitting region is illuminated continuously by a hidden intense beam, 
discussed by Rees et al. (2000), it can achieve the NEI (Tg < T^) state by the photoionization 
process. Even if this process realize, the mean energy of the hidden-beam photons which 
illuminate the line emitting region should not be largely greater than the edge energy of 
9.28 keV, since the observed values of T^ were low of kT^ = 0.8]'^q;2 keV at the rest frame for 
ASCA and kT^ > 1 keV for Chandra respectively. 

We pay attention to the case (ii) of a rapid adiabatic expansion of a highly ionized 
hot plasma. This sort of mechanism has already been investigated by Itoh & Masai (1989) 
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for a supernova which explodes in the circumstellar matter ejected during its progenitor's 
supergiant phase. They show that when the blast shock breaks out of the dense circumstellar 
matter into a low-density interstellar medium, a rarefaction wave propagates inward into the 
shocked hot plasma. Then the hot plasma expands adiabatically and loses its internal energy 
quickly. The electron temperature Tg decreases, but the ionization degree T^ temporally 
remains high, because the recombination time scale becomes much longer due to the low 
density. The mean temperature of the shocked matter, kT ~ 100 (fg/lO^ cm s~^)^ keV 
where Vg is the shock velocity, drops by about two orders of magnitude for the density 
contrast (ratio) of the dense matter to the ambient medium of ~ 10^ in their hydrodynamic 
calculations. Therefore, the plasma can achieve the NEI (Tg ~ 1 keV, T^ ~ 100 keV) state 
naturally, and emit the strong RRC with the weak Ka line of -FRRc/-^iine ~ 4. 

If these mechanisms work, we can estimate an emission measure for GRB 970828 using 
the observed photon flux of the RRC (^rrc), the distance D and the integrated emissivity 
(eRRc) shown in the figure 2 with the cosmical abundance, 

nV ~ io68( D )2/ ^RRC w gRRC .-1 .2) 

^3Gpc^ 4.7 X 10-5^4 X 10-16^ ^ ' 

Although these n and V are coupled with each other and highly depend upon the specific 
model, we may conclude that the density of line emitting region is considerably high. 

The RRC and/or the line of a large EW suggest a recombing (Tg < T^) condition, which 
can be realized by the photoionization or rarefaction. It should be noted that most of X-ray 
afterglow did not show both emission features; NEI is not always the case. 
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Fig. 1. — Simulated emissivities, convolved with the energy resolution of the ASCA-SIS for 
the cases; (a) T^ = 10 keV, Tg = 10 keV (equilibrium), (b) T^ = 1 keV, T^ = 10 keV, (c) 
Tz = 15 keV, Tg = 2 keV and (d) T^ = 100 keV, Tg = 1 keV. The solid lines represent 
emissivity of continuum and the dotted ones are for emission lines. The figure 1-c and 1-d 
are simulated for representing the cases of GRB 991216 and GRB970828 respectively only in 
view point of the observed ratios of -FrrcZ-^Uiic- We do not intend to reproduce the spectral 
shapes, which mostly consist of a non-thermal component. 
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Fig. 2. — Ratios of the integrated emissivity of the RRC structure to the Ka hues: 
(-^RRcZ-^iinc) for several ionization states in 3-dimensional axis. We calculate the ratio of 
the emissivity in each T^ = 5 keV and Tg = 0.5 keV steps. The straight line labeled "equi- 
librium" indicates the condition of T^. = T^. The recombination process enhances only in 
the Te < Tz region. Within our calculated range, the observed ratio of -FRRcZ-^unc > 3.3 for 
GRB 970828 can be attained in the range of more than Tz = 80 keV, limiting the Tg = 1 keV, 
which is observed in the RRC structure. 
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Table 1: Intensity of Iron Line and RRC 



GRB 


z 


-feline 


o 


-^RRC 


kT 


Name 




(photons cni~^s^^) 


(keV) 


(photons cm-^s-^) 


(keV) 


970508* 


0.835 


(3.0 ±2.0) xlO-5 


— 


— 


— 


970828^ 


0.958 


< 1.5 X 10^6 


— 


1.7+6;4 X 10-5 


0.8^^:^ 


991216tt 


1.020 


(3.2 ±0.8) X 10-5 


0.23 ±0.07 


3.8 ±2.0 


> 1.0 


000214^ 


— 


(9 ±3) X 10-6 


— 


— 


— 


990123* 


1.600 


< 3.3 X 10-6 


— 


— 


— 


990704* 


— 


< 4.7 X 10-6 


— 


— 


— 



* Piro et al. 1999 ; the line intensity was variable. 
"I" Yoshida et al. 2001 ; the RRC was temporal. 

« Piro et al. 2000 

^ Antonelli et al. 2000 

* Yonetoku et al. 2000 
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